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Abstract 

We show that the nonhnear optical response reflects sensitively the electronic 
structure of transition- metal surfaces and interfaces, d and s electrons may 
contribute rather differently to the second harmonic generation (SHG) signal. 
This results from the different sensitivity of d and s electrons to surface and 
symmetry changes. Consequently, SHG for noble metals shows a by far larger 
dependence on the polarization of the incoming light than for transition metals 
like Fe, Co, Ni, in particular at lower frequencies. The theoretical results are 
compared with recent measurements. We conclude that the SHG yield is in 
addition to the nonlinear magneto-optical Kerr-effect a sensitive fingerprint 
of the electronic structure at surfaces and interfaces. 
78.66.Bz,73.20.At,42.65.Ky 
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I. Introduction 

Recently, the nonlinear optical response of metallic surfaces has been studied inten- 
sively. The theoretical analysis uses typically the free-electron approximation or the jellium 
model 0,0. Clearly, this must be extended for studying the nonlinear optical response of 
transition metals and for analyzing d- vs. s-electron contributions to optical second har- 
monic generation (SHG). Here, we present such a theory by treating the electrons using 
a tight-binding like model which has been successful in analyzing the nonlinear magneto- 
optical Kerr-effect (NOLIMOKE) [^. As a result of the different localization of s- and d- 
wavefunctions and associated different sensitivity with respect to surfaces and interfaces 
and symmetry changes, we find that the wavelength and polarization dependence of the SH 
light on the incoming light is rather different for noble metals such as Cu, Ag, Au and for 
transition metals such as Ni, Co, Fe, and Cr, in particular at lower frequencies. 

Denoting by Ia{(3 — SH) the yield of /5-polarized second harmonic light resulting from de- 
polarized incoming light, then quite general theoretical arguments predict that for transition 
metals involving d-electron transitions 

I,{p-SH)^I,ip-SH), (1) 

while 

Ip{p-SH):^h{p-SH)^0, (2) 

if only s-p electronic transitions are occurring. Hence, for noble metals Ip{p — SH) ^ 
Is{p — SH) ^ for light frequencies u which do not permit d-electronic transitions, e. g. 
huj < [ep — Here, ep is the Fermi energy and e'^°'^ refers to the upper edge of 

the occupied d-band. However, for hu > {sp — e^°'^) noble metals will behave more and 
more like transition metals. 

This behavior results from the fact that for s (and p) electrons one gets for the dipole 
matrix elements controlling the nonlinear optical response (s | x | s) ~ and 
{s \ z \ s) 7^ 0, and furthermore {d \ z \ d) > {d \ x \ d) for transitions involving 
(i-electrons. The coordinates z and x refer to directions perpendicular and parallel to the 
surface, see Fig. 1. This behavior of the nonlinear optical response is in agreement with 
recent experimental observations To understand the effect of the matrix elements on 
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the nonlinear optical response in more detail, note that the nonlinear optical response is 
given by 

= x^^Je.Ei, (3) 

where Pj is the i-th component of the polarizability and where the susceptibility tensor xlji 
is determined within an electronic theory using response theory . Thus, the SHG yield is 
determined by 

J.(/3-5if)~|Exg;p (4) 

ijl 

neglecting for simplicity at the moment Fresnel and transmission coefficients. The nonlinear 
susceptibility is now of the form 

= E h I I J I k") {k" \l\k) F,,,,.{lu, 2u) , (5) 

kk'k" 

where {k \ i \ k') is the usual dipole-matrix element involving electronic transitions between 
s, p, and d states k and k'. Note, as a result of inversion symmetry breaking at surfaces and 
interfaces s — > {s,p), s — > d, and d — > d transitions occur The function Fi:k'k"{^^) 
collects all other factors (energy denominators and Fermi functions) and gives essentially 
the doubly convoluted density of states 0. Decomposing XiJ with respect to the electronic 
transitions one gets 

X^l = + xSjT + xi?''*, (6) 

with Xiji^'^^ involving only s, p — > d transitions, while Xifi^ involves only s — > s transitions 
and Xifi"^ involves only d — > d transitions. It can be shown that due to the dipole matrix 

(2) 

elements the various contributions to the tensor elements xlji of different magnitude, 
see table I. This will simplify the analysis. 

Using now for a simplified treatment of the wave functions the atomic basis set of ref . , 
we calculate the dipole matrix elements. Results are listed in table I. Thus, we get typically 
for transition metals Xifi'^ ^ xlfi'^ ^ ^-^id for noble metals xlfi ~ xlfi'^ ^ if hu < {ep — f^JT"^)- 
As a consequence, one estimates for noble metals with hco < {ep — ^^^^ due to 

{s \ X \ s) ^ 0: 
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Ip{p - SH) > Is{p - SH) , Isip -SH)^0 



(7) 



Furthermore, for transition metals we estimate: 

UP-SH)r^\Y.X^^]'\\ (8) 

ijl 

and thus 

Ip{p-SH) ^ 3h{p-SH), (9) 

since {d \ z \ d) ~ {d \ x \ d), see table I. 

This behavior of the second harmonic intensity — SH) is in accordance with re- 
cent experiments and shows already clearly that the SHG yield reflects sensitively the 
electronic structure of noble metals and transition metals at surfaces. Of course, this is of 
great interest not only for surface studies, but also for the investigation of thin films and 
interfaces. Note, for a more-detailed study of the polarization dependence of /^(p — SH), 
with a varying from p to s, Fresnel and transmission coefficients have to be included to 
determine the SHG yield. 

In section II we present a simplified theory for the polarization dependence of the nonlinear 
optical response which may serve as a basis also for further extensive studies. In section III 
we present results and compare with recent experiments. 

II. Theory 

That the SHG yield is a sensitive fingerprint of the electronic structure at surfaces and 
interfaces becomes obvious from the polarization dependence of the SHG light intensity 
— SH). Here, a refers to the polarization of the incoming light. In order to determine 
more accurately and in detail this intensity Ia{l3 — SH) of the /5-polarized SHG light at 
frequency 2uj, note that the refiected light at frequency 2uj is given as shown by Bohmer et 
al. [^Pl for symmetry without magnetization by 
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tl sin^ 
'^fsipts cos V9 sin 
y 2fctpts COS sin y 

Here, $ and denote the angles of polarization of the reflected frequency doubled and of 
the incident light (see Fig. 1). /c,s are Fresnel coefficients and ts,p are the linear transmission 
coefficients. The complex indices of refraction at frequencies u and 2uj are n = rii + iki 
and N = n2 + ik2- The projections of the incident wave vector k on the spatial coordinates 
inside the medium are fg = sinO/n and /c = y^l — fg. The corresponding quantities 
for the reflected SHG light are Fg = sinQ/N and = y^l — F^, where 9 and 6 denote 
the angle of incidence of the incoming light and the angle of reflection of the reflected SHG 
light, respectively. The linear transmission coefficients are given by 0,^,0 

2cos^ 2cos^ „ 2cose „ 2 cos 9 



tr. 



T 



^ ncos^ + /, ' cos^ + n/e' NcosB + F^ 

The corresponding amplitudes Ap and Ag in Eq. (10) are 

and 

Ag = 



Tg 



cos e + NF^ 



2ttTp 
cos 



2nTg 



cos 6 



(12) 



respectively. From Eq. (10) one gets directly the intensity /($, (f) of the reflected SHG light: 
/($, ip) = I £'(^'^)($, ip) p. Thus, one obtains for p-polarized SHG with vacuum permittivity 
e and vacuum permeability fiQ 



Up-SH) = 2 1^ I {2^-) n Eoiiv) ^ A,((F,xa2/J. + iV2F,(xEU^ + xSl/f ))t>osV 
+ N'Fgxfltlsm'^)\\ (13) 



and for s-polarized SHG 



Us - SH) = 2,/^ I {2t-) n Eoioo) H AsX^^l2U,t, cos ip simp \^ . (14) 
V A*o c 

Here, Ia{P—SH) denotes the intensity of /3-polarized reflected light at frequency 2a; resulting 
from a-polarized incoming light {a — p,s,(p — t^/^)- In particular, one gets for p-polarized 
incident and p-polarized reflected SHG light 

I,{p - SH) = 2 If I (2^) ri E,{u) n A,tl{F,x^l2Uh + N^F^ixfLfl + xflfs)) ? , 

(15) 

and for s-polarized incident and p-polarized reflected SHG light 

hip - SH) = 2 /1^ I (2^^) E,{u) n A,N^F,xfLtl ? , (16) 

and furthermore for "mixed" -polarized (45° polarization) incident and s-polarized SHG light 

W4(« - SH) = 2M I {^) n £;o(a;) \'\ A,x^l2U,ts T ■ (17) 

\l Ho c 

The response in the latter two conflgurations is determined by one single tensor element 
xfJx ^^<i xSL' respectively. Then we get for an incidence angle of 45 ° for the ratio 
Ip(p - SH)/Is{p - SH) the result 

I,{p - SH) _ , (FcXii2/J, + N^F,xf)j^sftl , 



(18) 



Is{p-SH) ' iV^F^xiftf 
where the contribution from x^zxx been neglected in the numerator. Note, this tensor 
element has been experimentally determined from /(^=,r/4(s ~ SHG) and has been found to 
be at least three times smaller than the other two tensor elements (see below). Eq. (17) 
holds also for noble metals. 

It follows now from Eq.(18) that s and d electrons contribute rather differently to the 
intensity of SHG light: 

First, for noble metals hke Cu, Ag, Au one gets using Eq. (16) due to 

) 



xfii ~ (s I i I s) (s I j I s) (s I Z I s) and (s \ x\ s) = 



Is{p-SH) = 0, (19) 

if the light frequency uj is such that hu < {sp — ^d^""^) implying that no transitions 
involving d-states occur. Note, x^^Jx Eq. (16) involves the matrix elements {s \ x \ s) and 
these are nearly zero. Hence, for noble metals 
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Up - SH) ^ ^'"^ 

if the optical response does not involve d state transitions. Thus, Eqs. (18) and (20) suggest 
a very sensitive test of the quantity [ep — e^°'^) and of many-body contributions to this. 
Comparison with results from bandstructure calculations may help to identify correlation 
effects. 

If the frequency u increases, then for hu > {ep — ^T'"'^) ^Iso s — > d and d — > d 
transitions contribute and 

Lip - SH) ~| j:tiFsX^' f + I EtlFsxfJ: f + .... (21) 

ijl ijl 

Hence, for increasing frequency u, Is{p — SH) is no longer (nearly) zero, increases and 
Is{p — SH) — y Ip{p — SH) as for transition metals. 

Secondly, for transition metals like Fe, Co, Ni, Cr, we estimate in contrast to noble 
metals for all frequencies u, that 

Ip{p-SH) ^ Is{p-SH). (22) 

For example, for the Ni(OOl) surface, using hu = 2.056 eV (A = 603 nm) and for the 
optical constants |0 rii = 1.98, ki = 3.92, n2 = 2.02, and k2 = 2.18, we find using Eq. (17) 
for 9 = 45° 

if we allow the relative phase of the tensor elements xf}z ^'^d xi^L to vary between and 2-k. 
This result gives already a reasonable account of the experimentally observed ratio of 1.25. 
Note, the estimate of the ratio Ip{p — SH)/Is{p — SH) was obtained only on the basis 

(2) 

of symmetry considerations for the matrix elements in xlji and by approximating these by 
atomic orbital matrix elements. Of course, for discussing systematically the nonlinear optical 
response of different noble and transition metals we have to include density of states effects 
resulting from the function Fkk'k"{'^). How much this matters can be seen from results for 
xfxxlx^llx. This ratio should be equal to one if Fkk'k"{^) plays no role. This corresponds 
to the case of perfect Kleinman symmetry [|ri|]. However, note that in xfji the indices 
j, I refer to uj photon transitions, while i belongs to the 2uj photon transition . This fact 



can be rather important, for example, for analyzing differences between noble metals and 
transition metals. For transition metals Xxji rather unfavorable for the SHG yield, since 
s — > s transitions at frequency 2u may occur and reduce the intensity I^^^/4^{s — SH) 
due to (s I a; I s) <^ {d \ x \ d), see for illustation Figs. 2c and 2d. In contrast, for noble 
metals XxJi is favorable for the SHG yield, if the transitions at frequency uj involve s — > s 
transitions only and if d — > s transitions at frequency 2uj come into play. The situation is 
illustrated in Figs. 2a and 2b. Thus, the deviation from Kleinman's symmetry of the ratio 
xixx/ X^xzx is interesting, a fingerprint of density of states effects, and reflects to some extent 
the contribution of d electrons to the optical transitions and the degree of localization of the 
Fermi-level electrons. Experimentally this ratio ranges from less than 0.013 for Al (with no 
d electrons) over 0.085 for Cu (where the upper d-band edge is 1.9 eV below the Fermi-level) 
and 0.33 for Au (where the upper d-band edge is 1.8 eV below the Fermi-level) to 0.35 for 
Ni (with a large density of d- electrons at the Fermi level) . For the determination of the 
ratio xfxx/x^L see remarks below. 

As can be seen from the large range of values for the ratio \^^^Zsh) i^ -^l' (^3), it is 
very important for a quantitative analysis to include not only the symmetry properties of 

(2) 

the dipole matrix elements, but also the phases, since the tensor elements Xiji are complex 
quantities. As a result, important interference contributions may occur and change the 

(2) 

results obtained from a coherent superposition of the partial intensities given by | xlji I 
alone. For a full analysis one has to perform an electronic calculation of Xiji as done 
previously by Pustogowa et al. 0. To compare theory and experiment, the latter has to 
be performed for several optical geometries to allow for a complete determination of the 

(2) 

absolute values and phases of xlji- 

In determining the phases at least approximately we may proceed as follows. To see 
how the two relative phases between the three tensor elements xizl? X^xL ^-nd xixx entering 
Ia{p — SH) come into play, we rewrite Eq. (13) as 




Xzzz 



with 



8 



y(2) 

^xzx I A.XZX I 1(^1 

X^ZZZ X.ZZZ 



and 



(2) v(2) 

e"^\ (25) 



^(2) - I ^(2) 

XZZZ X.ZZZ 

.(2) , , , A2) 



Now, the two absolute ratios | ^^jw- I and I ^fff I and the two relative phases and (x>2 remain 

zzz A.ZZZ 

to be determined. As can be seen from Eqs. (16) and (17), the intensities I^=T,/i{s — SHG) 
and Is{p — SHG) depend exclusively on one single nonlinear tensor element each, x^xzx 
X^zxx-i respectively. Thus, the intensity ratio 

I^=-K/i{s — SH) 
Is{p-SH) ' 

which may be taken from experiment, gives directly the absolute ratio | X^xzxlxfxx I- Using 
for example for /<^=7r/4(s — SH)/Is{p — SH) the experimental value and the appropriate 
linear optical constants of Ni one gets 

1^1= 1/0.35 ^ 2.8. (26) 

Xzxx 

Similarly, we calculate (1/0.013) for Al, (1/0.085) for Cu, and (1/0.33) for Au. 

To obtain the remaining relative phases </?i and </72, and the absolute ratio | X^xlxlx^zzz I 
one may use the results for the p-polarized SHG-yield for incident p polarized light. Note, 
according to Eq. (15), Ip{p — SH) involves a superposition of all three independent nonva- 
nishing tensor elements xfLi Xxzx xiaL- Keeping Ffefe/fe//(a;, 2a;) fixed, we may calculate 
all surface transition matrix elements entering xfji within an basis set of atomic 4s, 4p, 
and 2>d wave functions. Using these matrix elements (see table I), the decomposition of the 
nonlinear tensor xfji in Eq. (6) and Eq. (8) together with the partial p polarized SHG yields 
for transition metals in table II, gives us for Ni the absolute ratio 

Y^^) 3 

I /^xzx I " 



Xzzz 

It follows then from Eqs. (16) and (17) that the remaining absolute ratio is 



(27) 



y(2) y(2) (2) 3 

•^zxx I I A.ZXX II A.XZX I n OCT n 01 ^'OQ^ 

-(27 1 - 1 -(2r II -(2r I - ■ 5 - (28) 
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Using then Eqs. (24), (27) and (28) and the experimental results for Ia{p — SH) one can 
determine the remaining two relative phases of the three tensor elements. We find good 
agreement with experiment, if we use (for Ni) 

^(2) ^ (2) il.957r 

and 

^(2) ^ (2) i0.5057r (oq) 

A.ZXX Azzz^ ' I 

see results presented in Fig. 3. Thus, we have completely determined the absolute ratios 
and the relative phases of all non- vanishing elements of the nonlinear tensor x-^/ for surfaces 
of fee crystals. 

It is important to note, that this procedure shows how the polarization dependence 
of the SHG yield can be used to determine the various tensor elements of the nonlinear 
susceptibility x-^/. Furthermore, it shows also how in turn the microscopic susceptibility 
Xiji determines the polarization dependence of the SHG yield. 

Another interesting fingerprint of the electronic structure is the behavior of I^{s — SH). 
While it follows straightforwardly from symmetry and matrix elements, see Eq. (14), that 

Ip{s-SH) ^ Is{s-SH) ^ 0, (30) 

the behavior of the ratio 1,^=450(5 — SH)/Is{p — SH), in particular the observation 

fIip='i5°(,S - SH) ^ /y=45°(-S - SH) 

ICu > { TT TTfT^ — )Au [oi-) 



Isip-SH) ^ Isip-SH) 



(2) 

requires a more careful theoretical analysis, including in the discussion of xlji ^ot only the 

(2) 

dipole matrix elements, but again also the energy denominators. Note, xlji can be written 



as 



xg? = const. E {k I ^ I k'){k' \ j \ k"){k" \ I \ k) ^^Ijl ^ , (32) 



k,k',k" 



where x^^'* denotes essentially the linear susceptibility, but without matrix elements 



Thus, 



^(1) _^(i) 

i^k,k',k"(^>2cu) = const. ^^'^ 
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Hence, Ia{s — SH) reveals in addition to the matrix elements the behavior of the function 
Fkk'k"{^,'^^) involving the density of states. Note, the u dependence must follow from 
Fkk'k"{^, 2^;) and will be more general than the uj dependence deduced from the jellium. 

To understand now the result given in Eq. (31), note that (1^=450(5 — SH))cu > 
(1^=45° (s — SH))au, since the density of p states above ep is larger for Cu, while {Is{p — 
SH))cu < {I sip ~ S H)) Aui since the density of d states at Tiui below e^? is larger for Au [p!2| . 

Note, the ratio 

/y=45°(g - SH) 

Is{p-SH) 

probes essentially the density of d electrons at an energy of hu below the Fermi-level and 
of p states above ep, see for illustration Fig. 2. As a consequence, it follows from Eqs. (16) 
and (17) that the ratio xfxx/x'^L which describes the possibility of d electron excitations by 
optical transitions increases from the jellium type metal Al via the noble metals Au and Cu 
to the transition metal Ni [§,0. 

Thus, it is interesting to note that the combination of the symmetry properties of the 
(real) transition matrix elements and the density of states effects will account for a quan- 
titative understanding of the SHG yield. Furthermore, since the tensor elements xlji ^ire 
complex quantities, their phases cause important interference contributions to the coherent 
superposition of the partial SHG yields resulting from the single tensor elements alone. An 

electronic theory along the lines of Pustogowa et al. is required for the calculation of the 

(2) 

full complex SHG tensor xl/i to explain the experiments quantitatively and in detail. It 
is clear from Eqs. (10) - (17) that experiments have to be carried out for several optical 
geometries to sort out for comparison with theory the tensor elements and to allow for a 
complete determination of their absolute values and phases. 

The physics described by our theory is evident also from the results shown in table I. 

(2) 

Here, we have estimated the dipole matrix elements which determine largely Xiji by using 
atomic s{p) and rf-electron wavefunctions. Decomposing x^^j according to Eq. (6), it follows 

Xffj » X^, (33) 

since {s \ x \ s) ~ 0, and 
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(2,d) > (2,d) 

A,zzz A,zxx 1 V"^/ 

since {d \ z \ d) ~ {d \ x \ d). Generally, for transition metals, xfjf'' is largest by far. For 
noble metals, xfjf^ begins to dominate over xfji"^ huo increases and becomes larger than 
the energy distance of the d states from the Fermi energy Sp, hu > {ep — e]^""^). Further- 
more, approximately Ip{p — SH) is essentially given by and x^xL^ whereas Is{p — SH) is 
determined by xlz*!^ alone. In the case of s-polarized SH light /^(s — 5//) = Is{p—SH) = 0. 
And I^p=■,^/4:{s — SH) is given by x^xL alone. Hence, using table I, we find quite generally 
that: 

(a) for noble metals such as Cu, Au, Ag 

Ip{p-SH) > I,{p-SH), ifhu < {ef - ^r") 
Up-SH) ^ I xgr^ ? + I xSr^ ? , if /^u; > {sf - eT^). 

lp{p - SH) > Jtrans.met.(^p _ gjj^ ^ y^^^^^ 



(b) for transition metals such as Ni, Co, Fe, Cr 
I{p-SH) ^\x^f \\ 

and 

Ip{p-SH) ^ 2h{p~SH). 
This completes then the theoretical analysis of the polarization dependence of the SHG yield 
and shows how it reveals the electronic structure at surfaces and interfaces. 



III. Results and Discussion 

To demonstrate the different polarization dependence of the SHG yield for noble metals 
and transition metals, we show in Fig. 3 numerical results for the polarization dependence 
of the p polarized SHG yield I^pij) — SH). Using Eq. (24), we calculate I^p{p — SH) for 
Ni and Cu where the p-polarized SHG- yield at p-polarized incident polarization has been 
normalized to unity. The angle characterizes the polarization of the incident light. We 
show results for Ni and Cu at angle of incidence ^ = 45 ° corresponding to the experimental 
situation The absolute ratios and relative phases of the tensor elements are listed in 

table HI as well as the used linear optical constants at frequencies uo and 2uj. 

First, we refer to results for Ni given by the topmost curve. These are typical for 
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transition metals. For Ni we used an incident light wavelength of 603 nm. The results are 
in perfect agreement with experiment Despite the fact that we have fitted the phases 

ifi and to experiment this is remarkable, since the ratio | xi^x/xE^ I has been obtained 
using only the matrix elements calculated with atomic wave functions. As mentioned the 
line shape is characteristical for transition metals. Similar results are expected for Cr and 
for Fe and Co. The SHG yield for s and p polarized incident light is nearly the same in 
transition metals due to the isotropy of the transition matrix elements for d electrons 

{d \ z \ d) ~ {d \ X \ d). 

Slight intensity minima occur at polarization angles of 45° and 135 °. These arise from the 

(2) 

interference of the tensor elements xlji the Fresnel factors, see Eq. (13). The intensity 
profile is always symmetric around = 90° and for transition metals approximately 
symmetric around (f = 45° and (f = 135°. 

Secondly, to demonstrate the typical frequency dependence of the SHG yield as a function 
of incident polarization for noble metals (N. M.), results for Cu are given in Fig. 3 by the 
curves (a), (b), (c) . Note, the incident light polarization dependence of the SHG yield is 
found to be very different for noble metals such as Cu than for transition metals (T. M.). 
The (p = 90 ° symmetry survives. The p-polarized SHG yield falls off to one fourth or less 
at incident s-polarization compared to its (normalized) value at incident p polarization, due 
to the strong anisotropy of the transition matrix elements for s electrons 

{s \ z \ s) ^ {s \ X \ s) ~ 0. 

The Cu curves (a) and (c) refer to wavelenghts exciting and non-exciting the Cu d electrons, 
respectively. The results are in very good agreement with experimental ones for the frequency 
dependent SHG yield obtained for Cu by Petrocelli et al. [Q. While the curve (a), which is 
for an incident light wavelength of 516.6 nm, still shows shallow minima at 63 ° and 117 °, 
the curve (c), which is for an incident light wavelength of 652.5 nm, falls off monotonously 
to zero at incident s-polarization (90 °). These minima are due to the phases of the complex 
nonlinear tensor elements giving rise to interference contributions as well as due to the 
complex optical constants. These interferences will be absent if the d band of Cu can 
no longer be excited for larger wavelenghts. Note, the shape of curve (a) would also be 
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characteristic of Au at 603 nm as is observed in the experiment by Bohmer et al. 
These results indicate that 



{h{p-SH))M.M. ^ {Is{p-SH))t. 



M. 



for increasing frequency uj as claimed previously on theoretical grounds. 

Note, that experimentally Is{p — SH) is larger for Au than for Cu. The experiments 
are performed for hu = 2.056 eV. Then, hu > {ep — according to bandstructure 

calculations by Papaconstantopoulos |jl2| yielding {ep — e™"^) ~ 1.8 eV for Au and 1.9 



eV for Cu. Hence, one would expect for both metals Cu and Au already larger values for 
Is{p — SH). However, note that due to many-body band narrowing [ep — £^3*"^) is for both 
metals somewhat larger than what results from bandstructure calculations. However, since 
{ep — e^^'^^Au < i^F — £d^"'^)cu one also expects in agreement with experiment 



mp-SH))Au > (Isip-SH)) 



Cu 



at energy huj = 2.056eV. 

The curve (b) corresponds also to Cu at a wavelength of 652.5 nm as curve (c). For 
this curve (b), however, only the optical constants have been adapted to this wavelength, 
whereas the nonlinear tensor elements of wavelength 516.6 nm were kept as in curve (a). 
At (f = 90°, this curve reaches the same height as the Cu curve (a), but it shows deep 
minima at <y9 = 54° and at (f = 126°. Thus, this curve (b) clearly demonstrates that 
for the correct theoretical description of the frequency and polarization dependence of the 
SHG yield the frequency dependence of the nonlinear tensor elements is as essential as the 
frequency dependence of the optical constants at frequencies u and 2lj. 

In summary, the polarization dependence of the SHG-yield shows characteristic curve 
shapes which allow for a clear distinction between the electronic structure of noble metal 
and transition metal surfaces. As claimed, SHG is simply a very sensitive probe of the 
electronic structure. 

One very interesting application of our theory for the frequency dependence of the SHG- 



yield is to transition metal oxides |]T3|. Note, late 3d transition metal oxides, in particular 
NiO, are charge transfer insulators. NiO has a gap of 4 eV separating the filled oxygen 2p 
band from the empty upper Hubbard portion of the Ni 3d band, see Fig. 4 for illustration. 
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These transition metal oxides are of particular interest due to the strong electron correlations. 
Thus, the frequency dependent SHG could give some important information about the energy 
gap of the transition metal oxides and the width of the upper Hubbard band and thus about 
correlations. Fig. 4 illustrates the electronic structure of NiO, for example. 

Accordingly, one expects for example for NiO a.t hu < A, where A refers to the charge 
transfer gap, no SH yield. As a consequence, oxidization of the Ni surface (in air) would not 
essentially change the experimental results for Ni, since the Ni below the oxidized surface 
would mainly contribute to the SHG yield. Similarly, this conclusion can also be drawn for 
Cu. Actually, for Cu, the irrelevance of the oxide layer was shown by Petrocelli et al. 0] 
from the agreement of the experimental data in air and in ultrahigh vacuum. Furthermore, 
since the SHG-yield would disappear iifiu} ~ A and fiu} > W, nonlinear optics could be 
used to determine the width W of the upper Hubbard band. 

We have neglected hybridization between the oxygen 2p states and the majority 3d states 
of Ni. However note that at ep and above only minority d states are present. 

In view of the different magnetic properties of the NiO-vacuum interface and of the 
Ni-NiO interface, the nonlinear magneto-optical Kerr-effect (NOLIMOKE) may provide es- 
sential additional information regarding the electronic structure and the location of the 
electronic contribution to the SHG signal. Very interestingly, the nonlinear magneto-optical 
Kerr effect should vary characteristically with light frequency u, since the ci-electron spin po- 
larization varies throughout the c?-electron bands of Ni, whereas no nonlinear Kerr-signal is 
expected for antiferromagnetic NiO. Thus, in the range of optical frequencies, only minority 
d-electrons of Ni at the Ni-NiO interface contribute to the NOLIMOKE signal. 

Finally, it is of interest to apply our theory also to surface and interface states having 
generally different s,p and d characters and spin-polarization. 

We have neglected contributions to the SHG yield of higher than electric dipole order. 
However, we estimate them to be smaller by a factor of ~ (1/137)^ = 5 ■ 10~^ for the 
magnetic dipole and of {2nao/X)^ ^ 2.5 • 10^'' for the electric quadrupole contribution. 
Here, a denotes the fine structure constant, and is Bohr's radius. The resulting smallness 
cannot be compensated by the skin depth factor. The implied interface sensitivity of SHG 
is supported by the agreement of our theory with experiment. 
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The uj dependence of the SHG yield follows from the function Fkk'k"{^, 2cij), see Eq. (32). 
It could be interesting to sort out this in limiting cases in order to compare with results 
from continuum electrodynamics and jellium approximation 0]. 

In summary, our theory shows that the nonlinear optical response at transition and noble 
metal surfaces and interfaces is a sensitive fingerprint of the electronic structure. Both, 
matrix elements and the density of states affect characteristically the intensity of the SHG 
2uj reflected light. While a more detailed numerical analysis is desirable, it will not change 
our qualitative arguments based on general physical facts. 
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TABLES 

TABLE I. Non-dipole transition matrix elements at the surface involved in the calcu- 

(2) 

lation of xh/ and I{SH), respectively. 



mat rix (^Ic^mont 


A'aluo 


(4s 1 X ± 1 4s) 





(4s z 1 4s) 


as 
Z 


(4pm=±l 1 Z 1 4:Pm=±l) 


69 as 
16 Z 


(4pm=0 1 Z 1 4pm=o) 


69 as 
8 Z 


{Sdm=±2 1 Z 1 3dm=±2) 


105 as 
64 Z 


(3(i„i=±i z 3dm=±i) 


105 OS 
32 Z 


{3dm=o z 3dm=o) 


105 as 
32 Z 


{Sdm=o -2 4s) 


0.7005^ 


(4pm=±i \x±iy\ 4pm=o) 


69 as 
8%^ ^ 


{^Pm=0 \x±iy \ ■ipm=Tl) 


69 as 
8V2 ^ 


{3dm=±2 \ X ± iy 1 3drn=±l) 


105 as 
32 Z 


{3dm=^i \x±iy\ 3dm=^2) 


105 as 
32 Z 


{Sdm=±2 \x±iy\4:s) 


0.7005 /l5as 
8 y 2 Z 


{Sdm=±i \x±iy\ 3dm=o) 





{Mm=o \x±iy\ 3dm=^i) 
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TABLE II. Values for the ^^polarized partial SHG yields I'ff''''^^ - SH) of d and s 
electrons and interference contributions. 



partial SHG yield 


value 


I^^l^ip-SHG) 


25(^)6(^)6 


lf}'£{p-SHG) 


9(^)6(^)6 


ll'}t\p-SHG) 




ifl'^iv-SHG) 





I?lt%-SHG) 


0.015152(^)6 




0.052908(^)^ 



TABLE III. Values for the absolute values and relative phases of the nonlinear tensor 
elements xlji sind for the optical constants rii, ki of Ni and Cu corresponding to the four 
curves of Fig. 3 (A denotes the wavelength of the incident light). The optical constants 
for Ni are taken from ref. [10] and for Cu from ref. [14]. 





Ni (A=603 nm) 


Cu (A=516.6 nm) 

curve (a) 


Cu (A=652.5 nm) 
curve (b) 


Cu (A=652.5 nm) 
curve (c) 


1 Xxzx \ 
1 (2) 1 

y ?:?:?: 


0.60 


0.363 


0.363 


0.0363 


1 Kzxx \ 
\ (2) 1 
"Xxzx 


0.35 


0.085 


0.085 


0.085 


phase ifi 


1.945 TT 


1.02 TT 


1.02 TT 


1.02 TT 


phase (p2 


0.505 TT 


1.34 TT 


1.34 TT 


1.34 TT 


ni{u) 


1.98 


1.12 


0.214 


0.214 


ki{uj) 


3.92 


2.60 


3.67 


3.67 


n2{2oj) 


2.02 


1.53 


1.34 


1.34 


k2{2io) 


2.18 


1.71 


1.81 


1.81 
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FIGURES 

FIG. 1. The polarization and geometry of the incoming uj and reflected 2a; hght, respectively. 

and $ are the polarizations of the incident light and the reflected frequency-doubled photons. 

= 0° corresponds to p polarization and ip = 90° to s polarization. denotes the angle of 
incidence. The crystal axes x and y are in the crystal-surface plane whereas z is parallel to the 
surface normal. 

FIG. 2. Illustration of second harmonic generation from noble metals [(a) and (b)] and tran- 
sition metals [(c) and (d)]. For noble metals, in case (a) no d electrons can be optically excited, 
which in contrast is possible in case (b). In case (c) for transition metals, predominantly d electrons 
contribute to the SHG yield, whereas in case (d) the excitation starts from the s band. 

FIG. 3. Numerical results for the SHG yield Iip{p — SH) for Ni and for Cu, curves (a), (b), and 
(c), using Eq. (24). ip refers to the polarization of the incident light. For input parameters see 
text and table III. 

FIG. 4. Scheme of the electronic structure of NiO and Ni. W is the width of the upper Hubbard 
band, A is the charge transfer gap, U denotes the on-site Coulomb repulsion, eg and tg are the two 
and threedimensional represenstations of d electron states in cubic environment. 
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